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Numerical  study  of  interaction  of  vortex  structures  in  plasmas  and  fluids' 

V.Yu.  Belashov,  R.M.  Singatulin 

Kazan  State  Power  Engineering  University,  5 1 Krasnosel'skaya,  Kazan  420066  Russia 

The  results  of  numerical  study  of  evolution  and  interaction  of  the  vortex  structures  in  the  continuum , and, 
specifically , in  plasmas  and  fluids  in  the  2D  approach , when  the  Euler-type  equations  are  valid , are  pre- 
sented. The  set  of  the  model  equations  eidtxl  = dyM/B,  e^d^i  = -dxM/B,  3?p  + v • Vp  = 0 , 

v = -(zx  V\|f)  / B , Ayj  = ~p  describing  the  a continuum  or  quasi-particles  with  Coulomb  interaction  mod- 
els, where  p is  a vorticity  or  charge  density  and  t jf  is  a stream  function  or  potential  for  inviscid  fluid  and 
guiding- centre  plasma,  respectively,  and  H is  a Hamiltonian,  was  considered.  For  numerical  simulation 
the  CD  method  specially  modified  was  used.  In  terms  of  vortex  motion  of  fluids  the  results  of  numerical 
experiments,  specifically,  showed  that  for  some  conditions  the  interaction  of  vortexes  in  continuum  may  be 
nontrivial  and,  as  for  the  “classic”  FAVRs,  lead  to  formation  of  complex  forms  of  vorticity  regions,  for 
example,  the  vorticity  filaments  and  sheets,  and  also  can  ended  to  formation  of  the  turbulent  field.  The  un- 
dertaken approach  may  be  effective  in  studying  of  the  atmospheric  and  Aljven  vortex  dynamics,  and  also 
useful  for  the  interpretation  of  effects  associated  with  turbulent  processes  in  fluids  and  plasmas. 


L Basic  equations 

In  this  paper  we  study  numerically  the  evolution  and 
interaction  of  the  vortex  structures  (so-called  FAVRs 
[1])  in  the  continuum,  and,  specifically,  in  plasmas  and 
fluids  in  2D  approach,  when  the  Euler-type  equations 
are  valid.  In  general  case  the  set  of  the  model  equations 
describing  the  a continuum  (inviscid  incompressible 
fluid)  or  quasi-particles  (charged  filaments  aligned  with 
a uniform  field  B)  with  Coulomb  interaction  models  is 
the  following: 

eidJxi  = 3y.H/B,  *idtyi  =-dxH/B,  3V  = 9/3v, 
3,p  + v-Vp  = 0,  v = -(z  xV\}f)/B,  (1) 

Ay-/=-p 

where  et  is  the  charge  per  unit  length  of  the  filaments 
or  the  strength  (circulation)  of  discrete  vortex,  <J>  is  a z- 
component  of  vorticity  C,  or  charge  density  p,  and  \jf  is  a 
stream  function  or  potential  for  2D  flow  of  inviscid 
fluid  and  guiding-centre  plasma,  respectively,  and  H is  a 
Hamiltonian.  Note,  that  in  the  continuum  (fluid)  model 
£=1  in  the  Hamiltonian  eqs.  (1).  Function  /=  0 for  the 
continuum  or  quasi-particles  (filaments)  with  Coulomb 

interaction  models  [2],  and  f = k2\\r  for  a screened 
Coulomb  interaction  model  [3].  We  will  consider  here 
only  a case/=0,  and  generalization  of  our  approach  for 

/ = k2y  is  rather  trivially. 

For  numerical  simulation  the  contour  dynamics  (CD) 
method  [1],  to  some  extent  modified,  was  used.  This 
method  gives  a possibility  not  only  to  observe  evolution 
of  single  contour,  but  also  to  study  the  interaction  be- 
tween vortexes  having  different  symmetry  order  (differ- 
ent modes).  Let  us  demonstrate  it  in  the  next  paragraph. 


2.  Numerical  results  and  discussion 

Let  us  consider  some  results  of  numerical  simulation  in 
terms  of  the  vortex  motion  of  the  inviscid  incompressi- 
ble fluid,  as  more  visual.  In  general,  to  study  the  evolu- 
tion of  vortex  structures  with  different  symmetry  orders 
it  is  necessary  to  insert  a small  amplitude  disturbance 
r = R0  [1  + e cos  (ma  - £5m  t)]  (where  R0  is  a conditional 
radius,  £ is  an  eccentricity,  m is  symmetry  order  (mode), 
a is  an  angle  and  8>m  = (m  - 1)/  2 ) to  the  constant  vor- 
ticity circle  region.  But,  accounting  that  the  results  of  evo- 
lution for  one  and  two  vortices  with  different  m were  de- 
scribed in  detail  in  [4],  let  us  stay  on  results  on  interac- 
tion of  vortices  and  consider  the  most  simple  case  of  cir- 
cle vortices  when  m=  1 and,  therefore,  = 0 . For  two 
vortices  the  result  of  the  interaction  depends  on  sign  of 
vorticity  £ ("polarity")  and  the  distance  5 between 
boundaries  of  vortices.  For  two  vortices  having  opposite 
polarities  for  initial  distance  8 = d where  d is  a vortex 
diameter  we  observed  that  the  vortices,  rotating  in  op- 
posite directions,  move  in  the  same  direction  and,  prac- 
tically don't  interact  independently  on  value  of  8.  For 
the  vortices  having  opposite  polarity  the  result  of  evolu- 
tion depends  essentially  on  8.  So,  for  rather  small  8 the 
vortices,  on  a level  with  rotation  about  their  own  axes 
and  around  of  their  common  center,  interact  forming  a 
common  vortex  region  which  consists  of  the  vorticities 
of  more  small  scales.  For  rather  big  8 the  vortices  on  a 
level  with  rotation  about  their  own  axes  rotate  around  of 
common  center,  at  this,  their  interaction  is  reduced  to  a 
cyclic  change  of  their  shape  (so-called  "quasi-return" 
phenomenon  [1]  is  observed).  In  our  numerical  experi- 
ments we  have  found  that  critical  initial  distance  divid- 
ing these  two  types  of  interaction  §cr  = 3d  1 4 . 
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To  study  the  interaction  of  the  linearly  disposed  vortices 
in  more  details  we  have  consider  the  problem  for  four 
vortices  being  at  initial  time  along  one  line.  The  simula- 
tion of  the  interaction  at  linear  disposition  of  four  vor- 
tices one  can  observe  that  for  rather  big  and  equal  initial 
distance  between  vortices  the  evolution  leads  to  forma- 
tion of  two  vorticity  regions  as  a result  of  more  strong 
interaction  of  the  "outer"  vortices  with  the  "inner"  ones. 
At  this,  the  interaction  of  forming  pairs  is  similar  to  that 
of  two  vortices  case.  In  case  8,  = d I 2 we  observed  the 
formation  of  a complex  vortex  structure  which  consists  of 
many  vorticities  of  more  small  scales.  The  further  evo- 
lution of  this  structure  leads  to  formation  of  complex 
turbulent  field.  Let  us  note  that  in  last  case  we  can  also 
see  that  the  interaction  between  outer  vortices  is  more 
strong.  This  fact  can  be  explained  by  the  fact  of  more 
strong  "attraction"  of  outer  vortices  to  the  "center  of 
mass"  of  the  vortex  system  because  the  outer  vortex  is 
attracted  to  the  center  by  three  other  vortices,  and  the 
inner  vortex  is  attracted  to  the  center  by  two  vortices 
and,  in  opposite  side,  by  one  outer  vortex.  To  test  this 
statement  in  the  next  series  of  numerical  experiments 
we  have  arranged  outer  and  inner  vortices  on  different 
initial  distances.  As  a result,  we  observed  the  formation 
of  vortex  structure  from  two  inner  vortices. 

In  the  next  series  of  experiments  we  studied  the  interaction 
between  the  vortices  disposed  at  initial  time  in  the  comers  of 
appropriate  equilateral  figures.  At  this,  the  following 
results  were  obtained.  In  case  of  evolution  of  three  vor- 
tices with  opposite  signs  of  £ being  at  initial  time  in  the 
comers  of  triangle,  we  have  obtained  that  a pair  of  them, 
having  opposite  polarities,  behaves  as  well  as  pair  of 
vortices  with  opposite  polarities  in  two  vortices  case, 
and  third  vortex  does  not  participate  in  interaction  al- 
most, practically  independently  on  value  of 
5;  (/  = 1,2,3).  The  similar  character  of  interaction  is 

observed  for  four  vortices  with  opposite  signs  of  ^ being 
at  t - 0 in  the  comers  of  square. 

The  character  of  interaction  of  three  and  four  vortices 
having  the  same  polarities  depends  essentially  on  the 
distances  between  them  like  in  the  two  vortices  case. 
For  example,  for  5 = d / 2 < 8cr  we  observed  that  three 
vortices  are  rotated  forming  one  big  vortex  which  con- 
sists of  many  vorticities  of  more  small  scales.  Similar 
picture  is  observed  for  four  vortices  being  at  / = 0 in  the 
comers  of  square  on  distances  smaller  critical  one  from 
another. 

3.  Conclusion 

So,  we  have  presented,  as  more  visual,  some  results  of 
numerical  simulation  of  eqs.  (I)  in  terms  of  the  vortex 
motion  of  the  inviscid  incompressible  fluid.  But,  as  we 
noted  in  the  beginning,  the  set  of  eqs.  (1 ) with  / *=0  may 
describe  also  the  quasi-particles  with  Coulomb  interac- 
tion model.  At  this,  the  results  presented  above  can  be 
easily  extended  on  the  2D  simple  system  where  the 
plasma  is  represented  by  charged  filaments,  aligned 
with  a uniform  field  B,  that  move  with  the  guiding- 


centre  velocity  ExB ! B2 . Moreover,  the  approach  un- 
dertaken can  be  useful  and  also  for  other  2D  continuum 
models  when  /*0  in  the  Poisson  eq.  (1).  They  can 
describe  the  vortices  or  filaments  with  the  non-Coulomb 
interaction.  In  the  last  case  it  is  assumed  that  ions  move 
with  the  guiding-centre  velocity  but  electrons  have  a 
Boltzman  distribution,  at  this,  the  additional  term 

/ = k2\\f  will  describe  the  Debye  screening  [3].  An- 
other plasma  model  that  can  be  investigated  using  the 
described  approach  is  the  Hasegawa-Mima  model  [5], 
its  equations  can  also  be  put  in  form  (1)  by  way  of  in- 
clusion of  ion  polarization  drift  through  the  ion  equation 
of  motion  [3]  dt  v = (el  A/)(-V<p  + vxB) . As  to  a hy- 
drodynamic fluid  model  corresponding  to  a screened 
interaction  that  such  model  having  form  (1)  has  been 
proposed  in  [6]  to  describe  the  Earth's  atmosphere  with 
the  equation  of  motion  for  the  horizontal  atmospheric 
flow  dfx  = -g¥h  + R\xi)  where  h is  the  atmospheric 
depth  and  R is  the  Coriolis  force.  Another  possible  ap- 
plication can  take  a place  in  the  study  of  the  problems 
associated  with  the  dynamics  of  the  Alfven  vortices  in 
the  ionospheric  and  magnetospheric  plasma  [7]. 

In  conclusion,  let  us  note  that  the  applications  of  the 
dynamics  of  different  types’  vortex  structures  will  not 
rest  the  problems  of  study  of  the  processes  discussed 
above,  in  the  power  units.  Such  investigations  are  now 
one  of  the  most  perspective  trends  in  some  technical 
fields,  for  example  in  the  problems  of  study  of  the  vor- 
tex motions  in  working  chambers  of  the  varied  types’ 
power  units,  where  the  working  substance  may  be  wa- 
ter, vapor-gas  mixtures  and  plasma  [8].  At  this,  the  most 
actual  problem  is  the  investigation  of  the  vortex  motions 
in  the  spatial  regions  where  there  are  the  flows  with 
velocity  shear  (for  example,  near  the  walls  of  the  boiler 
of  thermal  power  units  and  the  chambers  of  the  gas  tur- 
bines), and  also  the  study  of  the  vortex  motions  in  mag- 
netized plasma  without  dissipation  (controlled  thermo- 
nuclear fusion  systems)  where  the  electric  currents  can 
lead  to  formation  of  very  complicated  forms  of  the  vor- 
tex regions. 
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